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Abstract 

This paper concerns an application of the fuzzy 
structures analysis (FSA) procedures of Soize to 
prototypical aerospace panels in MSC/NASTRAN, 
a large commercial finite element program. A 
brief introduction to the FSA procedures is first 
provided. The implementation of the FSA methods 
is then disclosed, and the method is validated 
by comparison to published results for the forced 
vibrations of a fuzzy beam. The results of the 
new implementation show excellent agreement to 
the benchmark results. The ongoing effort at 
NASA Langley and Penn State to apply these 
fuzzy structures analysis procedures to real aircraft 
panels is then described. 

Introduction 

In 1986 Christian Soize was at the Office Na- 
tional d 'Etudes et de Recherches Aerospatiales 
(ONERA) 1 * in Chatillon, France and introduced the 
concept of fuzzy structures 2,3 ’ 4 ’ 5 ' 6 in the dynamic 
analysis of large scale structures. Having nothing 
to do with the use of fuzzy sets or fuzzy logic, 
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a fuzzy structure is composed of (1) a well char- 
acterized master structure and (2) attached fuzzy 
substructures whose details are either unknown or 
are imprecisely known. In the context of a tradi- 
tional finite element model, the master structure 
would be that part of the structure amenable to 
conventional FEs. The fuzzy substructure would 
be those parts of the structure inaccessible to tra- 
ditional modeling: missing degrees of freedom, fine 
scale structural details, attached equipment, etc. 
An implementation of the fuzzy structures concepts 
in FEs as originally conceived by Soize is called 
Fuzzy Structures Analysis (FSA), a term coined by 
Ruckman. 7 

Although there have been some attempts to 
implement alternative fuzzy structure approaches 
in the context of analytical modeling, 8,9,10,11,12 the 
original FSA methodology of Soize has great attrac- 
tiveness since it can be performed using existing FE 
models: no new special elements need be generated. 
What has been lacking, however, is an implementa- 
tion of FSA in a mainstream FE code. 

This paper describes an original implementa- 
tion of Soize’s FSA in the well known commercial 
code MSC/NASTRAN, 13 specifically its solution 
sequence 108, the dynamic frequency response mod- 
ule. Using NASTRAN’s Direct Matrix Abstraction 
Program (DMAP) environment and a MATLAB 14 
module to provide appropriate input parameters, 
the new implementation allows one to include fuzzy 
substructures along with any structural model al- 
ready existing in NASTRAN. The new DMAPs 
have been benchmarked against published results 
for vibrations of beams with attached fuzzy sub- 
structures, and these benchmarks are reported here 
for the first time. This paper also describes an on- 
going application of the new FSA implementation 
to two representative aircraft fuselage panels under 
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study at the NASA Langley Research Center. 

This paper is arranged with the following sec- 
tions: Firstly, the background of Soize’s fuzzy 

structures analysis is provided, focusing on both 
the overall concepts of the method and the required 
mathematical formulation. Secondly, this paper 
describes the specific implementation of the fuzzy 
structures analysis procedures in MSC/NASTRAN 
and MATLAB. Next a benchmark solution is de- 
scribed that was used to test the procedures. Two 
aircraft panels are then described which are the 
focus of our ongoing study in applying the fuzzy 
structures analysis procedures to real aircraft pan- 
els. 

Note that the acronym FSA is used for fuzzy 
structures analysis throughout this paper. Also 
for consistency with Soize’s works, the e+ tujt time 
convention is used. 

Background of FSA 

Motivation 

Soize introduced fuzzy structures for modeling 
systems with a huge number of degrees of freedom 
with a reduced size model. A good example of a 
fuzzy structure, a system too complicated for con- 
ventional modeling, would be the electronics under 
the cockpit of a commercial airliner. This mechan- 
ical system cannot be modeled directly because it 
is so complex. However, it must be modeled some- 
how since it contributes to the overall dynamics and 
structural acoustic response of the airplane. This 
particular fuzzy structure is important for interior 
cabin noise upon landing since it may be directly 
above or near the forward landing gear of the air- 
liner, and the landing gear vibrations can travel 
along the fuselage and cause acoustic radiation into 
the cabin. 

Another aim for FSA procedures is to account 
for uncertainty of small scales in structural dy- 
namics models. One source of uncertainty is in the 
variability in the construction of large airliners. An- 
other source of uncertainty is in the finite element 
modeling procedure itself, since no individual can 
possibly explicitly incorporate into a deterministic 
finite element model all of the fine structure in a re- 
alistic airliner model. One idea put forth by Jeffrey 
Cipolla 15 in 1995 was to think of FSA as a large 
eddy simulation (LES) for large dynamics prob- 
lems in analogy with how fluid dynamicists handle 
fluid scales too small for explicit representation in 


a direct numerical simulation (DNS). 

The computational cost is low for FSA com- 
pared with the conventional Monte Carlo modeling 
of uncertainties in finite element simulations. In 
1986 Soize showed one could achieve similar results 
between a FSA and a traditional Monte Carlo model 
of an uncertain structure but with many fewer FE 
solves. 2 * This computational savings, along with 
keeping only a minimum number of DOFs in the 
FE model, is the great advantage of FSA. 

Previous work 

The algorithm for implementing Soize’s FSA 
for complicated structures has been available in the 
literature for many years. A somewhat complicated 
algorithm to be described briefly below, the proce- 
dure involves making a few runs of a conventional 
structural dynamics finite element analysis program 
with a large number of carefully constructed right 
hand sides (forcing functions). The results obtained 
using each of the forcing functions is then handled 
using standard statistical techniques. 

In 1994 the Car derock Division of the Naval 
Surface Warfare Center tasked one engineer, Chris 
Ruckman, to travel to France to investigate Soize’s 
fuzzy structures methods in detail. One of Ruck- 
man’s goals was to determine the feasibility of 
implementing the fuzzy structures procedures in 
NASTRAN, and he wrote a short report 16 on how 
he thought one could accomplish this. He also 
published an even shortened version of the report 
as a paper at conference. 7 After the work was 
written up, changing priorities and reduced bud- 
gets at NSWC and the Office of Naval Research 
stopped further development of implementing the 
fuzzy structures procedures in NASTRAN. Ruck- 
man left NSWC, and he moved on to other areas. 
However, Ruckman correctly predicted that a fuzzy 
structures analysis in NASTRAN was possible. 

At about the same time the company now 
known as BBN Technologies of Cambridge, MA 
also was tasked by the Office of Naval Research 
to include a fuzzy finite element capability in their 
program SARA, a coupled fluid and structure finite 
element solver. 17 This development was accom- 
plished in 1994 using the Pierce/Sparrow/Russell 8 
and Strasberg/Feit 12 fuzzy structure models instead 
of Soize’s model. Because of reduced budgets, how- 
ever, no detailed testing or validation of the BBN 
implementation, to the best of our knowledge, was 
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ever made. SARA is not used in the commer- 
cial aircraft industry. However, since NASTRAN 
is more widely in use, the present authors decided 
that the prudent course was to implement Soize’s 
procedures in NASTRAN, putting them into the 
mainstream for testing and validation by NASA 
Langley, industry, and university researchers. 

Philosophy 

The general philosophy of FSA is to model 
as much of structure/fluid with conventional finite 
elements as possible. This portion of the structure, 
accessible to conventional modeling techniques is 
called the master structure. To add the influence 
of additional fine detail or uncertain attachments 
to the master structure, one modifies the boundary 
conditions of the master structure. See Fig. 1, where 
the Z' s indicate boundary impedances affecting 
the master structure. This addition is called the 
fuzzy substructure. By only incorporating the 
details of the fine or uncertain structure through 
the boundary conditions only, the total number 
of degrees of freedom (DOFs) in the model is 
unchanged. This helps keep the size of the finite 
element model from increasing to a point where 
conventional FE solvers are unusable. 

It should be noted that the fuzzy substructure, 
in the present anatysis, only affects the master 
structure through the FE nodes. One only places 
the fuzzy substructure at nodes where one expects 
the uncertainty or fine structure to exist. 

Fuzzy impedances 

Each additional boundary condition in the 
fuzzy substructure is modeled as the impedance 
of a one degree of freedom attachment. See Fig. 2 
where the mass fi is connected to the master struc- 
ture through a stiffness K and a damping constant 
C and where u) p = (K//i)^ is the natural frequency. 
Here £ = n is the critical damping constant as- 
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sumed to be small, (0 < £ < 1). Assuming the 
support has a normal velocity x = v = ve lut with 
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for the impedance of the attachment. 

To make the oscillator fuzzy/uncertain the 
following assumptions are made about the mass, 


damping, and modal density n of the oscillator; 

M = + U) 

f = £(i + y a ) (2) 

n = n(l + Yj) 

where _ indicates mean value, and where the Y t 
are independent random variables, in practice small 
fractions of unity. 18 

A one DOF oscillator having modal density 
may be unfamiliar, but it simply is equivalent to the 
oscillator having uncertainty in natural frequency 

Incorporation into the dynamic model 

In the following paragraphs U indicates a col- 
umn vector of generalized displacements, F a col- 
umn vector of generalized forces, and square brack- 
ets [ ] indicate matrices. A conventional structural 
dynamics vibration problem would be formulated 
for a particular angular frequency u) using the FE 
equation 


ioj ([^master structure]) ^ ^ (^) • (^) 

Adding the effect of the fuzzy substructure consists 
of modifying the dynamical equation at appropriate 
nodes as 

ii*J ( [^master structure] T [^fuzzy]) F F (td) - (4) 

Decomposing the fuzzy substructure into its mean 
and random components one has [Z( uzzy ] = [Z] + 
[Z ran d] yielding 

ibJ ([Zdct] 4* [Z ran d]) U — F (a;) (5) 

for the dynamical system, where [Zdct] = 
[Zmastcr structure] T [Z]. Soize theil pulls OUt [Zdctji 

i(jj [Z det ] ([/] — [T]) U = 7 (w) , (6) 

an equation which can be solved by a perturbation 
analysis. Here [T] = -[Z dc t] _1 [Z ran d] has entries 
assumed small compared to unity. Soize defines 


M 

[Zrand] = X>[Zt] (7) 

1=1 

where the [Z{] are the contributions to [Z ran d] 
by each random aspect of the fuzzy. There is 
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only 1 non-zero element in each [Z,]. There are 
M = L x J x d contributions to the sum: 

L = number of attachments 
J = number of directions (3 in general) 
d = number of fuzzy parameters (usually 3-4) 
and the Xi are independent, centered, uniform 
random variables. Note that the construction of 
the [Zi] y too lengthy to give here, is completely 
described in the works of Soize 2 ’ 3,4,0 as are other 
details concerning the solution procedure. 

Formally, if it could be constructed, the solu- 
tion of Eq. (6) is 

Because ||T|| < 1 one can write 

oo 

(w-mr 1 = [/]+£ [T\ k . (9) 

fc=i 

So Eq. (8) becomes 

w + £m MZdc.ir'FM (io) 

k=l / 

where 

M 

[T]=-[Z dct ]- 1 Ex«[^] • (11) 

i= 1 

The solution is constructed by expanding out 
Eq. (10): 

M M M 

u - U m + -£ + ■■■ 

i=l «i = l t2 = l 

( 12 ) 

Each of the U are solved for easily, and each 
of these solves is completely deterministic. The Xi 
are the only random variables in Eq. (12). 

Solution sequence 

The usual solution sequence follows from solv- 
ing Eq. (10) with increasing orders of the index k 
in these steps: 

0. FE solve with master structure and mean fuzzy 
only. 

1. Another FE solve of master structure with 
specially constructed right hand sides (RHSs) 
using results of 1st solve (Step 0). 

2. If necessary, another FE solve of master struc- 
ture using specially constructed RHSs using 
results of 1st and 2nd solves (Steps 0 and 1). 


3. Postprocessing of results gives mean and vari- 
ance of response. 

Here are additional details: 

Ste p 0: k — 0. The initial step gives by solving 
iu,[Z dcl }U (0) = 7(w) . (13) 

This is equivalent to solving Eq. (10) omitting the 
Slim completely. Note that the solution 
of Eq. (13) is completely deterministic. 

Step 1: k — 1. Get u[^ for 1 < i < M by solving 

iw[2dot] U[ l) = -iu[Zi]U {0) . (14) 

This is equivalent to solving Eq. (10) keeping only 
the first term of the sum - Usually 

Eq. (14) is implemented using 1 solve with M 
RHSs for each frequency Each frequency must 
be solved separately since [Z<j c t] turns out to be 
frequency dependent. Further, this k = 1 solve, 
again completely deterministic, usually is enough 
to represent Eq. (12) with sufficient accuracy. 

Step 2: k = 2. If necessary one further solves for 
U^] 2 for 1 < z'l < M and 1 < i 2 < M by solving 

i^[Z det ]U^l = -iu; [Z 12 }V^ . (15) 

This is equivalent to solving Eq. (10) keeping the 
first two terms of the SUTn - Here Eq. (15) 

is implemented using 1 deterministic solve with 
M 2 RHSs. Since a solve with M 2 RHSs could be 
prohibitively expensive, one hopes that this step is 
not necessary 

Step 3. Postp r ocessing: 

Now that the u[ l \ and if necessary U ^] 2 , 
components of Eq. (12) have been found, one can 
find the statistics of the solution. Ruckman 7,16 
elucidated this procedure. 

Clearly the solution is in the form of gener- 
alized displacements and the expectation (average) 

of U is simply U ^ obtained in step 0 of the so- 
lution sequence. However, the expectation of any 
magnitude squared value is 

£{M 2 } = t- (0) +£h (1) | ( 16 ) 

1=1 
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for fc = 1 when r is any scalar value such as 
displacement, velocity, acceleration, pressure, etc. 
The variance is further given for A: = 1 as 

M n p M j 

1=1 1=1 

M M o 

^EEM'.'V’)! ■ 

t= 1 j= 1 

(17) 

Implementation 

Overview 

To implement the solution procedure just de- 
scribed it was decided to use two consecutive NAS- 
TRAN DMAP runs, including modifications to the 
SOLUTION 108 (direct frequency response) mod- 
ule. The NASTRAN database capabilities would be 
used to minimize matrix recalculations between the 
two runs. The first NASTRAN DMAP run saves all 
the necessary values through a simple execution of 
the usual 108 sequence, performs precalculations, 
and saves needed information to the NASTRAN 
database. The second NASTRAN DMAP run ac- 
tually executes the k = 0 and k = 1 steps in the 
solution sequence described above as well as per- 
forming the solution postprocessing. An extensive 
MATLAB script is executed after the first NAS- 
TRAN run to precalculate information for the [Zi] 
matrices, required for the second NASTRAN run. 

The first DMAP code consists of 123 lines of 
DMAP inserted into usual SOLUTION 108 (direct 
frequency response) bulk data file. Here is where 
the usual 108 sequence is used. It also saves the 
database and counts the number of fuzzy DOFs. 
The second DMAP code contains over 650 lines of 
DMAP. It uses the NASTRAN database saved in 
the first DMAP, but is essentially a stand alone 
analysis. 

The MATLAB code, “makefuzzy,” provides 
inputs for the second DMAP code since NASTRAN 
can't perform 2-D numerical integration, required 
for generating the [Zi] matrices. 

Details on FSA DMAP 2 

Usually NASTRAN represents the mass, 
damping, and stiffness matrices straightforwardly 
as 

[Zdet] = -w 2 [M] + UJ [D] + [I<] . (18) 


However, in FSA the damping and mass matrices 
become frequency dependent as 

iw [Zdet] = ~ W 2 ([M] + [ttfuz HD 

+ iw ([Z? (w)] + [/f uz (w)]) + [A (w)] 

(19) 

as they now include the real and imaginary con- 
tributions of the fuzzy substructure, [Pf uz (^)] and 
[/fuzM], respectively. 

Hence the FSA DMAP 2 must solve for one fre- 
quency at a time. Note that this is not NASTRAN's 
default procedure. Usually MSC/NASTRAN at- 
tempts to perform the direct frequency response 
for all frequencies in solution 108 during the same 
solver pass. 

To explicitly clarify what has to occur, a 
solve for the 1st frequency would take the form 

iu;[Zdct(^i)] — F(^i), then one solves 

iw [Z dot (w, )] Z7, (1) = -iw [Zi (W! )] Z7 (0) ( Wl ) (20) 

for each i. Similarly, the corresponding solve for the 

2nd frequency would be zcj^Z^q^ Ml u {0) = F(w 2 ). 
then one solves 

iw [Zdet (wa)] u[ l) = -iw [Zi (w 2 )l u {0) (w 2 ) 

again for each z, etc. 

User procedure 

When running a typical FSA, the following 
procedure is followed: 

« Run FSA DMAP 1 

- (use scr=no command line option) 

- Input: nodes designated as fuzzy in a SET 

- Output: .DBALL database file 

- Runtime: about same as SOLUTION 108 

• Run Matlab Code “makefuzzy” 

- Input: frequencies for analysis, 

fuzzy parameters 

- Output: file fuzzyforphasela2.txt, 

matrix element values 

- Runtime: 5 min on a DEC Alpha work- 
station for 100 frequencies. 

• Run FSA DMAP 2 

- Input: fuzzyforphasela2.txt, 

“FSA DMAP U .DBALL file 

- Output: mean squared displacement, etc 
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- Runtime: 2 x (SOLUTION 108 runtime) 
4- time for RHS calculations. 

The number of RHSs — 9 x number of fuzzy 
nodes. This is because in the present implementa- 
tion only displacement DOFs (3) are assumed to be 
fuzzy, and only 3 fuzzy variables are assumed for 
each DOF. The 3 fuzzy variables are mass, critical 
damping, and modal density given in Eq. (2). The 
amount of deviation for these variables from their 
mean values are specified by the user as dispersion 
parameters A!,A 2 ,A 3 in the MATLAB procedure 
“make fuzzy. “ 

Testing and validation 

The resulting DMAP/Matlab procedures have 
been successfully applied to a benchmark prob- 
lem published by Soize, a fuzzy simply supported 
beam. 5 

Model parameters 

The base beam (master structure) parameters 
were taken directly from Soize J s work 5 . The beam 
mass was 20,000 kg, and the attached oscillators 
have a 0.002 damping rate and a mass of 1.8 kg. 
The dispersion parameters w r ere matched to be 
Ai = A 2 = A 3 = 0.001. 

However, the beam modeled here did have 
some differences to Soize’s beam. Soize’s model had 
2500 FE nodes and was modeled using a 0.05 Hz 
frequency resolution. He used 1213 1-DOF fuzzy 
attachments in the y (“up”) direction, with n = 
0.035, 0.043 for mean modal density on two distinct 
patches of oscillators. In the present work only a 
200 FE node model was employed using a 1 Hz 
frequency resolution from 1 to 100 Hz. Here 300 
1-DOF fuzzy attachments were distributed over the 
two patches specified in Soize’s problem, but they 
were in all 3 directions (a:, y, and z) and they 
all used the same value, 0.035, for mean modal 
density. The 300 fuzzy attachments were located 
at 100 fuzzy nodes, corresponding to 900 RHSs in 
the k = 1 FE solution. In summary for the present 
FSA DMAP 2, 100 k — 1 individual FE solves w T ere 
made, 1 solve for each frequency, with each solve 
containing the 900 RHSs. 

The run times on a NASA Langley Structural 
Acoustics Branch DEC Alpha workstation were 

• 1 minute — FSA DMAP 1 

• 5 minutes — “makefuzzy" Mat lab file 

• 6 h, 7 minutes — FSA DMAP 2 


Clearly the codes are not yet optimized for 
speed, but these numbers provide some evidence of 
the substantial time required for the procedure. 

Results 

The plot shown in Fig. 3 clearly indicates the 
similarity of the results for the mean squared ac- 
celeration at a point on the beam published as Fig. 
7(a) of Soize. 5 Here the solid line gives the fuzzy 
structures result and the dashed line gives the re- 
sult for the master structure alone. (The additional 
thin solid line in Soize’s Fig. 7(a) is one realiza- 
tion of a Monte Carlo simulation for comparison.) 
Sparrow’s beam model has a coarser discretization 
(200 FE nodes) and frequency resolution compared 
to Soize’s model (2500 FE nodes), but otherwise 
the plots agree very well in magnitude. A close in- 
spection of the two results shows that there is a 
slight discrepancy in the location of the resonances, 
and this again is due to the different discretizations 
of the models. These results clearly indicate that 
these new fuzzy structures procedures in DMAP 
are working correctly and are ready for application 
to the more realistic panel structures of interest to 
the aerospace and transportation industries. 

Note also in Fig. 3 the greatly constricted 
fluctuations of the structure’s acceleration with the 
addition of the fuzzy substructure. In essence, 
the fuzzy substructure acts to damp the master 
structure vibrations. 

Prototypical aircraft panels 

The present research at Penn State and NASA 
Langley is now focused on applying the new FSA 
procedures to realistic, prototypical aircraft panels. 
The investigation is centered on modeling the junc- 
tions between the longerons and skin and frames 
and skin in such panels as fuzzy/uncertain due to 
the inherent impreciseness of such joints. Note that 
this is a departure from Soize’s original applica- 
tion of the FSA procedures for uncertain attached 
equipment in large dynamical systems. 

Single longeron panel 

The first panel consists of a .04 inch thick 
aluminum skin with a single longitudinal stringer 
dividing the panel skin into two equal bays of 37.75 
by 7.625 inches. A single line of thirty 6-32 x 
3/8 screws attach the longeron to the panel skin. 
The panel is mounted in a rigid frame to simulate 
clamped boundary conditions. It is referred to as 
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either the single longeron panel, or the two-bay 
panel. 

Flat fuselage panel 

This second panel is representative of current 
aircraft construction but was manufactured without 
curvature to simplify the experimental and analyti- 
cal modeling. It is constructed of a 47- by 72-inch 
aluminum panel with a 0.050- inch skin with six 
equally spaced longerons and four equally spaced 
frame stiffeners. Single lines of rivets attach the 
stringers and frames to the skin. A bay is defined 
as a section of the panel skin that is bounded by 
the stringers and frames. The ba} r responses are 
the focus of the fuselage panel correlation efforts. 
For both the experimental and analytical studies 
on this panel, free-free boundary conditions were 
simulated. 

Progress in modeling 

Conventional FE modeling of the fuselage 
panel 19,20 has focused on different stiffener to skin 
attachment models and their effects on the pre- 
dicted bay motions. To characterize the panel 
dynamic response up to 1000 Hz, several methods 
of modeling the panel were examined. First, the 
required finite element mesh density of the panel 
skin was evaluated by performing a normal mode 
analysis of a single bay with clamped boundary 
conditions. The skin was modeled with linear plate 
elements. A mesh of 30 by 16 elements was found 
to provide a one-percent convergence on frequency 
and adequate spatial resolution to define the mode 
shapes through 1000 Hz. This resulted in 11682 
linear plate elements for the overall panel skin. 

For the FSA study, the stiffeners were modeled 
using one- dimensional beam elements with the ef- 
fective cross-sectional properties (area, inertia, and 
torsional constant) of the stiffeners. Offsets from 
the skin to the stiffener neutral axis are also in- 
cluded. The beam elements are created along a 
line consistent with the rivet line that attaches the 
skin to the stiffeners. The beam stiffener model 
of the panel contains approximately 60,000 degrees 
of freedom (DOFs). Similarly, the two-bay panel 
was modeled using a combination of linear plate 
elements for the skin and beam elements for the 
longeron. The resulting model of the two-bay panel 
has approximately 7900 DOFs. 

The FSA modeling effort is currently underway. 
Preliminary results indicate that one must be very 


careful in picking appropriate values of fuzzy mass, 
fuzzy damping, and fuzzy modal density as input 
parameters for the FSA procedures. If these studies 
are successful, plans are underway to apply the 
FSA procedures to more complicated prototypical 
structures including NASA Langley’s Aluminum 
Testbed Cylinder (ATC). 21 

Conclusions 

This paper has described an implementation of 
Soize’s fuzzy structures analysis (FSA) procedures 
in MSC/NASTRAN and MATLAB. The proce- 
dures have been validated by comparison to the 
results for a simply supported fuzzy beam already 
in the literature. Previous studies have shown that 
with FSA one can obtain similar results to a Monte 
Carlo simulation of uncertain attachments with a 
greatly reduced computational expense. The work 
now underway is aimed at the fuzzy modeling of 
the skin attachments in two prototypical aircraft 
panels, and detailed results are expected soon. 
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Figures 


Spray uncertainty onto inside of master structure . , . 




Figure 1: A FE model is made fuzzy by attaching appropriate boundary conditions to nodes. Each 

boundar}' condition is modeled as the impedance Z of an attached 1-DOF oscillator. Note the the total 
number of DOFs in the model is unchanged when the fuzzy substructure is added. 



x(t) 


Figure 2: The 1-DOF model used for each fuzzy DOF of the fuzzy substructure. The moving support, or 

master structure, moves with displacement x(f) and velocity £(£) = u(£) = ve lu>i . The resulting force due to 
the fuzzy 1-DOF attachment is F(t) = Fe liJt . 
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-40 


master structure (dashed), with fuzzy (solid) 



Figure 3: Mean squared acceleration at a point on Soize’s simply supported beam versus frequency. The 

dashed line shows the result for the master structure alone, while the solid line shows the result including 
the effects of the fuzzy substructure. The horizontal axis is the frequency in hertz and the vertical axis is the 
expected mean squared acceleration in dB relative to 1 m/s 2 . A I IIz frequency resolution is used. These 
results compare favorably to Fig. 7(a) of Soize’s work. 5 
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